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Abstract 

Some exact interactions between vibrational modes in systems with discrete symmetry can 
be described by the theory of the bushes of nonlinear normal modes (NNMs) [G.M. Chechin, 
V.P. Sakhnenko. Physica D 117, 43 (1998)]. Each bush represents a dynamical object conserving 
the energy of the initial excitation. Existence of bushes of NNMs is ensured by some group- 
theoretical selection rules. 

In [G.M. Chechin, et al. Int. J. Non-Linear Mech. 38, 1451 (2003)], existence and stability of the 
bushes of vibrational modes in the simple octahedral model of mass points interacting via Lennard- 
Jones potential were investigated. In the present paper, we study these dynamical objects by the 
density functional theory in SFq molecule which possesses the same symmetry and structure. We 
have fully confirmed the results previously obtained in the framework of the group theoretical 
approach and have found some new properties of the bushes of NNMs. 
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I. INTRODUCTION 


We consider nonlinear vibrations in Hamiltonian systems. Conventional, or linear normal 


Q. 


modes (LNMs) are exact solutions of dynamical equations in the harmonic approximation 
This approximation means, that potential energy of the considered system is decomposed 
into multi-dimensional Taylor series, and only quadratic terms are taken into account for 
obtaining Newton classical equations. LNMs cease to be exact solutions when we involve 
not only quadratic, but also some anharmonic terms. If the latter terms are sufficiently 
small, then it is possible to speak about some interactions between LNMs. In many physical 
problems, such interactions and the role of different anharmonic terms are studied. Taking 
into account the smallness of these terms, we can construct some approximate solutions for 
the considered nonlinear problem. As an example let us refer to studying the various types 


0 - 


of the phonon-phonon interactions in crystal physics 

In this connection one can ask, are there exist some exact solutions of nonlinear dynamical 
equations beyond harmonic approximation? 

Let we have a nonlinear system whose potential energy contains besides quadratic terms 
(therefore, the system admits the harmonic approximation) some anharmonic terms preceded 
by a common factor 7 . This parameter determines the strength of nonlinearity of the system. 


Lyapunov proved |3] that each LNM can be continued in 7 to obtain an exact periodic 
solution of the nonlinear system. He also gave a procedure to construct this exact solution 
up to arbitrary degree of the parameter 7 . As a result, one can find the set of N Lyapunov 
nonlinear normal modes (NNMs), where N is the dimension of the system. Unfortunately, 
the procedure of Lyapunov NNMs constructing is very cumbersome, and we can obtain 
solutions only up to a small degree of the parameter 7 . Moreover, this procedure usually 
converges only for very small values 7 . These factors prevent to use the Lyapunov modes in 
many physical problems. 


0 ). 


Another notion of nonlinear normal modes was introduced by Rosenberg 4] (see also 
These modes can exist for large values of the parameter 7 and even for the essentially non¬ 
linear systems whose potential energy does not contain any quadratic terms (such systems 
don’t admit the harmonic approximation). The number of NNMs by Rosenberg, if they 
exist, can be less or even greater than the full number of the system’s degrees of freedom. 

Unfortunately, Rosenberg modes can exist only in very specific nonlinear systems. Let us 
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consider this point in more detail for N-degrees-of-freedom mechanical system. According 
to definition, in the dynamical regime described by a given Rosenberg mode, all degrees of 
freedom, xft), vibrate as follows: 


Xiit) = a,if(t), i = 1..N. 


( 1 ) 


Thus, all Xi{t) at any moment t are proportional to the same time-dependent function fit) 
( a.t are constant coefficients). Actually, this means a separation of space and time variable^. 

Substitution of the ansatz (JT]) into differential equations of the dynamical system leads 
to the set of (N — 1) algebraic equations for obtaining amplitudes ai and to one differen¬ 
tial equation which determines time-periodic function fit). The latter equation is called 
“governing”. 

Rosenberg found several classes of mechanical systems for which Eq. ([T]) is fulfilled. The 
most important class is formed by systems whose potential energy is a homogeneous function, 
of arbitrary degree, of all its arguments. 

In fl Q], we have shown that existence of Rosenberg nonlinear normal modes can be a 
direct consequence of a certain discrete symmetry of the considered physical system. We refer 
to such modes as symmetry-determined nonlinear normal modes (SD-NNMs). Hereafter only 
this type of Rosenberg modes is considered. 


In 


)|, 7], we have also introduced the concept of bushes of NNMs for dynamical systems 
described by arbitrary group G 0 of discrete symmetry. Every bush possesses a certain 
symmetry group G which is a subgroup of the group Go (G C Go) and represents an exact 
solution of nonlinear dynamical equations. As a consequence, the energy of initial excitation 
of a given bush turns out to be trapped in this dynamical object until it loses its stability 
(then the bush transforms into another bush of lower symmetry and larger dimension). 
Every bush whose dimension m is greater than unity {m > 1) describes a quasiperiodic 
dynamical regime which is determined by m governing differential equations. In this sense, 
one can consider the Rosenberg mode as one-dimensional bush. 

The specific group-theoretical methods for constructing bushes of NNMs were developed 
-8]. It is essential that because of symmetry-related causes only a finite number of the 


m 


bushes of each dimension can exist in the physical system with a given symmetry group Go- 


1 The LNMs also satisfy Eq. CD with fit) = acos(wf + </>o), where a, do and w are amplitudes, initial phase 
and frequency, respectively. 
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For example, it was found {q-II] that for monoatomic chain with even interparticle interac¬ 
tions and periodic boundary conditions (such as FPU-/3 chain or electrical chain, considered 


m 


12]) only five SD-NNMs can exist (depending on the number of particles in the chain). 


All bushes of low dimensions for many mechanical structures with various types of transla¬ 
tional and point symmetry were found in R 3, Q- In particular, in jl3j, all SD-NNMs 
(one-dimensional bushes) were found for all possible mechanical structures with any of 230 
space symmetry groups. 

Note that construction of bushes can be done without any information about the type of 
interparticle interactions in the given physical system. One must know only its symmetry 
group and the geometrical structure. 

Searching for the complete set of modes entering the bush with a given symmetry, we deal 
with its geometrical aspect. This problem can be solved with the aid of group-theoretical 
method only. On the other hand, we deal with dynamical aspect of the bush when study 
time-evolution of its modes and the bush stability. In this stage of bush studying, we must 
know specific interactions between particles of the physical system. 


In 


14], we have studied geometrical and dynamical properties of one-, two- and three- 


dimensional bushes in a simple octahedral mechanical system which is depicted in Fig. [0 
It represents a regular octahedron with equal point masses at the vertices and with another 
point mass in its center. It was found, that 18 bushes of vibrational modes can exist in the 
case when the central particle is immovable. Among these bushes there are ID, 2D and 
3D bushes with point symmetry groups Oh, D^ and CJ,,, respectively. These bushes were 


Q 


studied in M not only as geometrical objects, but also as dynamical objects supposing that 
interparticle interactions in the considered system are described by arbitrary pair potential, 
U(r ), in particular, by the Lennard-Jones potential 

A B 


U (r) = 


„12 


r 6 ' 


( 2 ) 


Naturally, one can ask: “Are there exist some physical systems whose nonlinear dynamics 
can be described by the above simple mechanical model?” The more general question can 
be formulated as follows: “Can the concept of bushes of nonlinear normal modes and the 
methods of their studying be valid for nonlinear dynamics of real physical systemsD 

First of all, let us note that there are some molecules, whose equilibrium state corresponds 
to the mechanical model in Fig. [0 As an example, we choose the molecule SF (> . The exper- 
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imental investigation of nonlinear vibrations of such molecules represents great difficulties 


15] devoted to study interactions between normal modes 


m 


(see, for example, review paper 
simple molecules). 

On the other hand, the ab initio calculations based on the density functional theory 
(DFT) [l6-18], proved to be very effective and rather correct for studying molecules and 
crystals. Indeed, DFT allows to determine the geometrical structures of many microscopic 
objects up to 1% accuracy 18]. Therefore, it is very interesting to study nonlinear dynamics 
in the framework of DFT methods. In our recent paper [l9], we successfully realized such a 
program for studying discrete breathers in graphane. We found there that DFT leads to car¬ 
dinally different results for vibrations with large amplitudes than the methods of molecular 
dynamics based on the concept of mass points interacting via a widely used phenomenolog¬ 
ical pair potential. Indeed, the mass-point models cannot describe the complicated process 


of the polarisation o: 
Brenner potential 


the atomic shell during essentially nonlinear vibration. In 


20 ]. the 


21| was used for studying discrete breathers and the non-monotonic con¬ 


nection v(A) between breather amplitude A and its frequency v was obtained (it shows two 
changes between soft and hard nonlinearity). In contrast, with the aid of ABINIT pack¬ 
age 22|, we have obtained the monotonic function v(A), which demonstrates that only soft 
nonlinearity takes place for gap breathers in graphane. 

In this paper, we study nonlinear dynamics of SF e by DFT methods realized in ABINIT 
code 22], The paper is organized as follows. In Sec. 2, the geometrical aspects of nonlinear 
vibrations of SF 6 corresponding to the bushes with point symmetry groups Oh, D^h, C\ v 
are discussed. Sec. 3 is devoted to ab initio calculations of the dynamics of the above 
bushes. In Sec. 4, we comment on the validity of the bushes of the nonlinear normal modes 
in real nonlinear systems (using SFq molecule as an example) and discuss further possible 
investigations of nonlinear dynamics of physical systems in the framework of DFT. 


II. GROUP-THEORETICAL ANALYSIS OF NONLINEAR VIBRATIONS IN SF 6 


Group-theoretical analysis o 


was already fulfilled in Ref. 


nonlinear vibrations of the structure presented in Fig. [I] 
14], Here we reproduce some results of that work which are 


necessary for the further discussion. 

In any vibrational regime, the nuclei of the atoms, consisting the SF e molecule, displace 
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from the equilibrium positions depicted in Fig. [T] and we can speak about a displacement 
pattern at any fixed moment t. For every NNM, as well as for LNM, this pattern possesses 
a certain point symmetry group. In the conventional phonon-spectrum analysis (this is the 
case of small vibrations), the above patterns are determined by eigenvectors of the matrix 
of force constants. 


There exists Wigner theorem 23J about group-theoretical classification of the linear nor¬ 
mal modes. According to this theorem, the modes are classified by irreducible representations 
(irreps) of the symmetry group Go of the system in equilibrium. In this way, we can introduce 
the basis $ = {4>j\j = 1..N} in the space of all possible atomic displacements, constituted 
by the complete set of basis vectors of the irreps entering into the mechanical representation 
of the considered system. Therefore, any vibrational regime X(t) = {xi(t), x 2 (t ),... ,x n (t)} 
in this system can be decomposed into the above basis with coefficients depending on time t: 


N 


X(t) = £>i«<fe = (C(t),#). (3) 

j= 1 

In this equation, each term Cj(t)(f>j can be considered as nonlinear normal mode according 
to the definition ©• Indeed, the vector multiplier (pj determines the displacement pattern 
of all atoms, i.e. the space structure of NNM, while Cj(t ) determines time-evolution of the 
mode. However, for brevity, we often use the term nonlinear normal mode (or vibrational 
mode) individually for <pj, as well as for Cj(t). 

The basis vectors <pj correspond to different irreps T n of the group G o and, therefore, the 
displacement vector X(t) in ([3D can be written as the sum of contributions associated with 
individual representations of the equilibrium symmetry group G q: 


x (<) = X(C»«,«>|r„]). (4) 

Here, $[T n ] is the set of basis vectors of the irrep T n . 

According to Wigner theorem, the small vibrations of the molecule associated with the 
different irreducible representations T n are independent from each other. It means that if one 
excite (by using the appropriate initial conditions when solving linear differential equations) 
a dynamical regime X(£) corresponding to a given irrep T n , this regime can never leads to 
excitation of the modes belonging to another irreps in the decomposition (J3J) . Therefore, one 
can ask: “What will happen if we consider large and, therefore, nonlinear vibrations of the 
molecule?” The theory of the bushes of nonlinear normal modes starts from this question. 
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f. jf| (see also js], 


The answer was given in Ref. |6j (see also [8], devoted to discussion of the bush theory), 
ft turns out that there exist certain selection rules for excitation transfer from one mode 
to another. These rules are originate from some group-theoretical restrictions which can be 
written as a certain system of linear algebraic equations j6|. In particular, one can deduce 
from this system that excitation from the mode with the given symmetry group G can 
transfer only to those nonlinear normal modes whose own symmetry is higher or equal to 
G. The above selection rules lead to possibility for existence of bushes of NNMs. 

Each bush represents a set of NNMs that conserves the energy of initial excitation until 
it loses stability because of the phenomenon similar to the parametric resonance with some 
modes outside a given bush. This phenomenon occurs if amplitudes of some bush modes 
attain sufficiently large values (see details in BQ). 

Every bush possesses its own symmetry that is determined by intersection of all symmetry 
groups of its modes. As was already mentioned, when the given bush loses stability, it 
transforms into another bush with lower symmetry and with higher dimension. 

Let us consider the simplest bushes for nonlinear vibrations of SFq molecule using some 
results obtained in [lJ. In the equilibrium state, depicted in Fig. Q] the molecule SFq 
possesses point symmetry group Go = O All vibrational modes for this molecule, classified 
by irreps of the group Oh, can be found in Table. 3 in [bij] [three translational modes must 
be excluded since the central atom (S) is supposed to be immovable]. 

For the present consideration, we need explicit forms of the displacement patterns of 
NNMs (f) 1 , 0 2 , 0 3 , corresponding to one-dimensional irrep T 1; two-dimensional irrep T 5 and 
three-dimensional irreps T 10 . These displacement patterns are given in Table [D 


TABLE I. Displacement patterns of the NNMs in SFq molecule 


Irrep 

NNM 

Pattern 

IT 

01 

^(0,0, —1| — 1,0, 0 0, —1,0 1,0,0|0,1,0|0, 0,1 ) 

r 5 

02 

^=(0,0,21-1,0,010,-1,011,0,010,1,010,0,-21) 

ITo 

03 

^U(0,0, -2|0,0,1|0,0,1|0,0,1|0,0,1|0,0, -2|) 


In this table, for each fluorine atom, according to the numbering in Fig. [U we point 
out three coordinates x, y, z, which determines displacement of the nucleus of this atom 
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from the equilibrium position. One can see from Table 1 that the molecule shape in the 
vibrational regime, corresponding to the mode 0 1 represents, at any moment t, the regular 
octahedron. Its size vibrates in time becoming larger or lesser in comparison with the 
octahedron corresponding to the equilibrium state. This NNM is called “breathing” mode. 
It represents one-dimensional bush with symmetry group Oh- The symmetry group of the 
breathing mode, Oh, is higher than that of each other vibrational mode and, therefore, 
according to the above mentioned selection rules, the excitation from this mode cannot 
transfer to another modes. Therefore, the breathing mode will vibrate for arbitrary long 
time without involving into the dynamical regime any other vibrational modes. In other 
words, if we decompose X(t) for the breathing mode in accordance with the Eq. (J3J) , we 
obtain that right-hand side of this equation is reduced to one term only 

X(t) = ci(t)0i = a(*)<£i. (5) 

Here we have renamed the time-dependent coefficient c\ (t) as a(t). 

This means that breathing mode represents an exact solution to the nonlinear differential 
equations describing dynamics of our mechanical system for any type of interatomic interac¬ 
tions. The explicit form of the differential equation for the time-dependent coefficient a(t) 
can be found if one substitute Eq. ([5]) into the full original system of 3 x 6 = 18 nonlinear 
equations (remember that the central atom is immovable) taking into account the specific 
type of interatomic forces. In this way, we reveal that the above system of 18 equations is 
reduced to only one differential equation for a(t). Let us emphasize that a(t) indeed depends 
on interatomic forces, while the displacement pattern does not depend on these forces and 
can be found with the aid of the group-theoretical method only (see Table [J) . 

Thus, the breathing mode represents one-dimensional bush with symmetry group Oh and 
we denote it by symbol B[Oh]- In [Q, one can find an explicit form of the differential 
equation for a(t) in the case of arbitrary pair interatomic potential U(r), as well as for 
Lennard-Jones potential. 

Another situation appears when we excite the mode <p 2 by displacing fluorine atoms 
according to its pattern given in Table [B The own point symmetry of this mode is described 
by the point group G 2 = D 4h that is a subgroup of the equilibrium state symmetry group 
Go = Oh {D 4 h C Oh)- Displacement pattern associated with <f> 2 represents octahedron with 
square base, formed by atoms 2, 3, 4, 5 (see Fig. [TJ) , and the atoms 1, 6 situated on Z axis 


and displaced by the same distance from their equilibrium positions from the octahedron 
base. In Fig. [2J, we depict the above displacements of all six F atoms by arrows. In contrast 
to the case of excitation of the breathing mode, in the present case, the mode <p 2 cannot 
exist independently of all other modes. Indeed, its excitation leads to excitation of the 
breathing mode (f) { whose symmetry group G\ = Oh is higher than the symmetry group 
G 2 = D^h.. This conclusion is a result of the group-theoretical analysis only. Thus, we 
obtain the two-dimensional bush B[D 4 h] with the point group G 2 = D 4 h- From the general 
decomposition ([3]), we obtain for this case: 

x(t) = c 2 (t)0 2 + Ci(i)0! = b(t)(f > 2 + a(t)(f> 1 . (6) 


Here we have renamed ci(t) and c 2 (f) by a{t) and b(t), respectively. 

Substitution of the anzats ( 6 ) into 18 original nonlinear equations confirms that all these 
equations are transformed into two differential equations with respect to the unknown coef¬ 
ficients a(t) and b(t). These two time-dependent functions fully describe the dynamics of the 


two-dimensional bush -BjLL/jJi]. It describes quasiperiodic vibrations with two basic frequen¬ 
cies Ui and 0 J 2 corresponding to the modes a(t)<p 1 and b(t)<f > 2 (because of nonlinearity of the 
considered system, in the Fourier spectrum of such vibrational regime, one can observe not 
only the frequencies 04 , cn 2 , but also their different integer linear combination). 

Let us emphasize once more that no other modes present in the decomposition (|3lh i.e. 
Eq. ([ 6 ]) represents an exact solution of the original nonlinear equations whose coefficients 
Ci(t) = a(t) and c 2 (t) = b(t ) can be obtain as solution of two governing differential equa¬ 
tions depending on the specific interatomic interactions. In |l4|, we presented the explicit 
expressions of these governing equations for the case of arbitrary pair potential U[r ) and, 
in particular, for the Lennard-Jones potential. 

Let us note that in the vibrational regime described by the bush B[D 4 h] the atom con¬ 
figuration represents the octahedron with square base and with atoms 1 and 6 displaced by 
the same distance in opposite directions (this is a consequence of the presence in the group 
D 4 h horizontal mirror plane determined by the base of octahedron). 

The mode </> 3 (see Table |T|) , being excited at the initial instant t = t 0l generates three- 
dimensional bush B[C 4v \ whose full symmetry coincides with the group G 3 = C 4v of the 
mode (p 3 . 


1 Note, that a[t) in Eq. (J6J) is not in any relation with a(t) from Eq. ©! 
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The modes whose symmetry determines the symmetry of the whole bush we call “root 
modes”, while other modes which turns out to be excited automatically as a result of the root 
mode excitation we call “secondary modes”. Let us note, that in general the situation can 
be more complicated. Indeed, even for the original group G 0 = Oh there exist bushes whose 
excitation is possible only by exciting simultaneously several modes with different symmetry 
groups. In this case the symmetry of the whole bush is determined by intersection of the 
above symmetry groups. All nontrivial bushes of different dimensions (up to dimension 
equal to 8 ) were presented in fid} ]. 

In contrast to the group G 2 = D^., the horizontal mirror plane is absent in the group 
G 3 = C 4 As a consequence, displacements of the fluorine atoms 1 and 6 need not to 
be equal in magnitude and, therefore, the octahedron corresponding to the bush B[C± V ] 
possesses different heights dropped from the atoms 1 and 6 to the octahedron base formed 
by the atoms 2, 3, 4, 5. 

Thus, we obtain the three-dimensional bush B[C± V \ to which the following decomposition 
of X(t) corresponds: 


X(t) = c 1 (t)cj) l + cz(t)<t> 2 + c 3 (t)<j> 3 = a(t )(/) 1 + b(t)</> 2 + c(t) 4 > 3 . 


(7) 


Here we have renamed the coefficients Ci(t), c 2 (t), c 3 (f) by a(t), b(t), c(t), respectively. 
Substitution of the anzats (JH) into the original 18 nonlinear differential equations leads to 
three governing equations (all other original equations turn out to be equivalent to these 
governing equations) whose explicit form, for the case of pair potential U (r), are given in Q. 

For all above-discussed bushes, the symmetry of the vibrational state conserve, as well as 
the complete set of the bush modes, while amplitudes of these modes [they are described by 
coefficients cfit) in Eqs. (J5] ) .(J6 D ,(JTJ)] evaluate in time. 

We have considered three simplest bushes B[Oh], B[D ih \, B[C± V \ in the SF e molecule 
in the framework of the group-theoretical approach. The existence of these bushes was 
confirmed by straightforward numerical experiments using Lennard-Jones and Morse po¬ 
tentials. However, there is a principle question: “Can these bushes exist in real physical 
systems?” Indeed, during nonlinear vibrations of real atoms their electron shells polarized 
and it is not obvious that we can describe the complex physical process of such polarization 
(characterized by many degrees of freedom!) by simple mechanical mass-point models with 
pair potential interactions. There is a necessity to verify validity of the concept of bushes of 
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nonlinear normal modes and our group-theoretical methods, aimed at constructing bushes, 
using real physical experiments. Unfortunately, such direct experiments, as to the best of 
our knowledge, cannot be fulfilled with the aid of the present experimental technique. How¬ 
ever, we can use ab initio simulations based on the density functional theory for checking 
and deepening the theory of the bushes of NNMs. 


III. STUDYING BUSHES OF NONLINEAR NORMAL MODES IN SF 6 MOLECULE 
WITH THE AID OF THE DENSITY FUNCTIONAL THEORY 


The DFT is based on the very non-trivial theorem proved by Kohn and Hohenberg in 


16], 


Indeed, they proved that many particles wave function tH r i, r 2 , r 3 , • • • , Uv)> depending on 
3 N arguments and satisfying the stationary Schrodinger equation for a quantum system 
with N electrons, can be exactly expressed via the function of total electron density, p(r), 
depending on only three arguments r = ( x,y,z ). The function p(r) provides a minimum to 
a certain functional, J[p(r)], whose explicit form is not known yet, but one can be sure that 
it exists in principle. Some approximations were developed for this functional which turn 
out to be effective and rather accurate. 

In our work, we have used ABINIT code developed in 22] which realizes computational 
methods of DFT. This software package provides many useful facilities for ab initio calcula¬ 
tions of the structure and different physical properties of molecules, crystals and nanoclus¬ 
ters. 

Studying nonlinear vibrations of SFq molecule, we have used Born-Oppenheimer approx¬ 
imation for separating movement of heavy nuclei and light electrons, as well as local density 
approximation (LDA) and pseudopotential approach to single out core and outer electrons 
of the atoms consisting SFq molecule. 

Quantum-mechanical equations by Kohn and Sham 17] are used in the ABINIT code for 
describing dynamics of electrons, while classical equations are solved for nuclei with forces 
generated by electron shells at any time-step. For solving Kohn-Sham equations the basis of 
plane waves was used with maximal energy determined by cutoff energy -Ecutoff = 40 Hartree. 

Obviously, one can apply more complicated approximations such as GGA, etc. However, 
this paper is not aimed at calculations of physical properties which can be verified by real 
experiments (at the present time, we don’t know such experiments). We aimed at conhr- 
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mation of the general theory of the bushes of nonlinear normal modes in the systems with 
discrete symmetry using the molecule SFq as a simple example. In fact, description of this 
molecule in the framework of the density functional theory with some pertinent approxima¬ 
tions represents a mathematical model sufficiently closed to the real physical system. 

Let us discuss some numerical results on verification of the bush theory. 


A. The one-dimensional bush B[Oh\ 

First of all, we have calculated equilibrium configuration of SF ,6 molecule and have found 
that the edges of the regular octahedron is a 0 = 2.9828 Bohr. This value seems to be rather 
accurate since experimental value is equal to 2.9554 Bohr. 

At the next step, we excite vibrations of fluorine atoms by displacing their nuclei by some 
values, according to the displacement pattern of the breathing mode (j) 1 which represents 
the bush B[Oh] in Table U (see also Fig. 2a). In other words, we choose the initial shape 
of the regular octahedron, which can be determined by the edge a, larger [a( 0 ) > a 0 ] or 
smaller [a(0) < do] than the equilibrium size clq. Then we release the system, i.e. permit 
the molecule to evolve freely. In Fig. |3l we present the edge a(t) as a function of time for 
some initial values a(0). From this figure one can see periodic vibrations which turns out 
to be stable, at least, up to the time for which we have observed this dynamical process. 
The decomposition ((3]) of the corresponding vector X(f) = [xi(t) , X2(t) ,... ,Xis(t)] does not 
reveal (up to the numerical accuracy) any contribution of other modes besides the initially 
excited breathing mode cj) 1 . This is in full agreement with the group-theoretical analysis in 
the previous section. 

In Fig. [4l we depict the function z/(A), where v is the frequency of the breathing mode 
and A is its amplitude. As one can see from Fig. [4] the soft nonlinearity realizes in the 
breathing mode of SFq molecule. 


B. Two-dimensional bush B[D 4 f,] 

This bush can be excited by assigning certain values to its both modes at t — 0: < 22 ( 0 ) = //, 
cti(0) = v with zero initial velocities [di(0) = 0,a 2 (0) = 0]. This bush can appear only if 
a 2 ( 0 ) 7 ^ 0 (otherwise the one-dimensional bush B[Oh] will be excited). Indeed, the symmetry 
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of 0 2 i s l° w er than that of <fi 1 and, therefore, namely the mode 0 2 turns out to be the root 
mode for the two-dimensional bush B\D±h\. 

The most interesting way for the excitation of the bush B[Dy j] is the case fi 7 ^ 0, v = 0. 
In this way, we excite only the root mode 0 2 anc l have to see the excitation of the secondary 
mode 4 > 1 , because of the “force interaction” js] with the mode 0 2 , while all other vibrational 
modes of SF 6 molecule must not be excited. 

In Fig. [5J we present results of such excitation of the bush B[D^]. From this figure one 
can see that the amplitude £ii(f) = a{t ) of the mode 0 1; being zero at t — 0 , is gradually 
involved into the vibrational process. Then, with time evolution, we see how ai(f) acts on 
0 . 2 (f) and vice versa, or interaction between modes <f> 1 and 0 2 . 

Fig. [ 6 ] illustrates the dependence of max|oi(f)| on max|a 2 (t)|. Such dependence turns out 
to be important if we want to study the excitation of the secondary mode ai(f) 0 1 with the 
aid of the perturbation theory. 

Every bush can be treated as a closed Hamiltonian system. In the present case of the bush 
D Ahl we have such a system described by two dynamical variables ai(f) = a(f), a 2 (f) = b[t). 
In [71, 8 ], we study the “classes of dynamical universality” of the bushes. This concept 
appears naturally if one studies the decomposition of potential energy of different bushes 
into multidimensional Taylor series. Indeed, it turns out that many bushes in different 
physical systems, and even in the same system possess identical forms of the potential 
energy decomposition up to a certain degree. Because of this reason it is interesting to 
analyse decomposition of the potential energy of the bush B[D^f\. This can be done as 
follows. 

Let us choose a certain grid in the space of two variables /i and v. Each node of this grid 
determines the values of dynamical variables a(t) and b{t) of the bush B[D 4 0 at t — 0. The 
linear combination of the mode vectors 0 ,, 0 2 with the coefficients , corresponding 

to the chosen (i,j) node, determines a certain displacement pattern of the fluorine atoms 


Xjj( 0) = /iy0! + Vij(f> 2 (8) 

and, therefore, a definite initial configuration of the SFq molecule. 

On the other hand, using ABINIT code we can find the potential energy Uij of this 
configuration [for each node this program finds the electron density corresponding to the 
positions of the fluorine atoms determined by (| 8 j)]. Then one can construct the potential 
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energy U(a, 6 ) as a polynomial of a fixed degree in variables a, b: 


U(a,b) = Y,lmna m b n . (9) 

m,n 

We determine coefficients 7 mn by the least square method. As a result, the following potential 
energy U(a, b ) of the bush B[D 4h ] was obtained: 


U(a, b) = 0.134186 2 + 0.19269a 2 - 0.20382a6 2 + 0.042206 3 - 0.09507a 3 
- 0.04813a6 3 + 0.12396a 2 6 2 + 0.030856 4 + 0.03013a 4 . 


( 10 ) 


Here we point out only those terms of U (a, b) which are permitted by the symmetry group 
D 4 h, he. Eq. (fIU|) represents decomposition of U(a, b ) into the first polynomial invariants | 8 | 
of the group G = De¬ 
constructing U (a, b ) with the aid of the least square method, we have also calculated 
polynomial terms which are not permitted by the group De to verify that they are rather 
small (about 10 -3 — 10 ~ 4 ). 

Note that the similar expression for U(a, b ) we have already obtained in 14J applying the 
Taylor expansion of the exact potential energy for octahedral mechanical structure whose 
particles interact via Lennard-Jones potential. 

Now, we can write Newton equations for the two-dimensional system corresponding to 
the bush B[De\- 

a + 0.38538a = 0.28521a 2 + 0.203826 2 - 0.24791a6 2 - 0.12051a 3 + 0.048136 s , 

6 + 0.26835 b = 6(-0.126586 + 0.40763a - 0.123406 2 - 0.24791a 2 + 0.14440a6). 

From these equations, one can see “inequality of rights” of the root and secondary modes. 
Indeed, if we excite at t — 0 only secondary mode a(t ) by appropriate initial conditions 
[a( 0 ) 7 ^ 0 , a( 0 ) = 6 ( 0 ) = 6 ( 0 ) = 0 ], it will vibrate for arbitrary long time according to 
the equation a + 0.38538a = 0.28521a 2 — 0.12051a 3 , while the root mode remains to be 
unexcited: 6 (f) = 0. I 11 the contrary case, where only root mode 6 (f) is excited by the initial 
conditions [ 6 ( 0 ) 7 ^ 0 , 6 ( 0 ) = a( 0 ) = a( 0 ) = 0 ], the secondary mode, a(f), cannot be equal to 
zero, because there are “external” forces proportional to 6 2 (f) and 6 3 (t) in the r.h.s. of the 
first equation (fill) . 

n 

According to the bush theory (see [ 8 |) this fact is a straightforward consequence of the 
difference in symmetries of the modes a(f )</> 1 and b(t)<f> 2 which are Oh and D 4h , respectively 
(the group D 4h is the subgroup of the group Oh). 


( 11 ) 
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Solving Eqs. ([IT]) numerically with the aid of fourth order Runge-Kutta method we obtain 
some examples of the bush D 4 i, dynamics. Such example is presented in Fig. [TJ One ought 
to compare this figure with the similar Fig. ED which has been obtained by the ABINIT 
code. The comparison shows that Eqs. (HI]) , describing dynamics of the bush B[D 4h ] as 
an autonomous two-dimensional system, turn out to be rather correct for the considered 
vibrational amplitudes. 


C. Three-dimensional bush B[C 4v \ 

According to Eq. ([TJ) , the three-dimensional bush B[C 4v ] can be written in the form 

X(t)[Cy = a(t) 0 i + b(t)<f> 2 + c(t)<f> 3 . (12) 


Here, t/q, (f) 2 , 4 >3 are the vibrational modes from Table 1 which determine the special dis¬ 
placement patterns of the fluorine atoms of the molecule SF 6 , while a(t), b(t), c(t ) are 
time-dependent functions describing bush evolution in time. 

As was already discussed in Sec. EU the symmetry groups of the modes </q, <fi 2 , <p 3 are Oh, 
D 4h , C 4v , respectively. The following group-subgroup relation between these point groups 
exists: 

O h D D 4h D C 4v . (13) 


Thus, the mode <p 3 has the lowest symmetry group among three modes of the bush C 4v and, 
therefore, it turns out to be root mode of this bush. The point group of the whole bush 
B[C 4v ] coincides with that of the mode </> 3 . The modes <fi 9 and </q are secondary modes 
of the considered bush, and according to the bush theory [ 8 ] they must be automatically 
involved into the vibrational process as a consequence of the initial excitation of the root 
mode <f> 3 . 

To verify this prediction we have fulfilled simulation of SFq molecule vibrations with the 
aid of ABINIT code. Results of these simulations are represented in Fig.JSJfor three different 
values of the initial amplitudes c( 0 ) of the root mode <p 3 (two other initial amplitudes, as 
well as velocities of all three modes of the bush B[C 4v ] are assumed to be zero). This figure 
illustrates the excitation of the secondary modes <fi 2 and <p 3 because of their interactions 
with the root mode <p 3 . 
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Using the least square method, as it was described for the bush B[D^h], we have obtained 
the following polynomial expression for the potential energy U(a, 6, c ) of the bush B[C± V \ as 
the function of three variables a, b , c: 

U(a, 6 , c ) = 0.19280a 2 + 0.15057c 2 + 0.132546 2 + 0.211806c 2 - 0.23103ac 2 - 0.22387a6 2 
+ 0.043956 3 - 0.08708a 3 + 0.02771794903c 4 + 0.165226V + 0.031866 4 
- 0.2562a6c 2 - 0.05031a6 3 + 0.13336aV + 0.13581a 2 6 2 + 0.02 5 75a 4 . 

(14) 


Note that only polynomial invariants of the group G 3 = C^ v give large contributions to 
U(a , 6 , c) (all other polynomial terms obtained by direct least square method turn out to be 
rather small). 


It is very interesting to note that according to the theorems proved in {§], group-subgroup 
relation fjl3|l provides the special structure of the bush B[C± V \ dynamical equations. Indeed, 
the potential energy U(a, 6 , c) contains invariants which are linear in secondary modes, while 
they contain some degrees of root mode. We see such invariants in Eq. fjl4p : 


2 2 I 

c a, c 6 . 


(15) 


The Newton dynamical equations for the bush B[C 4 V ] as a closed Hamilton system can 
be written as follows: 

a + 0.3856a = 0.26125a 2 + 0.223876 2 + 0.23103c 2 - 0.10301a 3 

+ 0.050316 s - 0.27163a6 2 - 0.26672ac 2 + 0.2566c 2 , 

6 + 0.26 5 0 86 = -0.13 1 846 2 - 0.21180c 2 + 0.44775a6 - 0.127436 3 

(16) 

- 0.27163a 2 6 + 0.15093a6 2 + 0.256ac 2 - 0.330456c 2 , 
c + 0.30114c = c(0.46205a - 0.423616 - 0.11087c 2 - 0.26672a 2 
+ 0.512a6 — 0.330456 2 ). 

The invariants (fT5j) produce in the right-hand side of Eqs. (fTUj) forces which act from the 
root mode c(t) on the secondary modes a( 6 ) and 6 ( 6 ). Namely, these forces form the bush as 
unique dynamical object. The full set of its modes conserves in time, while their amplitudes 
a(t), 6 ( 6 ), c( 6 ) do change. 
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D. Nonlinearity in SFq dynamics 


In this paper, essentially nonlinear vibrations of SFq molecule are considered. Unfor¬ 
tunately, we do not know any experimental data that allow a direct comparison with our 
results. However, in the case of oscillations with small amplitudes these results should be 
consistent with those of the linear normal modes analysis. Let us discuss this point in more 
detail. 

Obtained in the harmonic approximation, linear normal modes represent sinusoidal os¬ 
cillations whose frequencies do not depend on their amplitudes, and all such modes are 
independent of each other. This means that if only one LNM was excited in the considered 
system, then it will exist for arbitrary long time without excitation of any other normal 
mode. 

Results of our ab initio simulations demonstrate a violation of these properties for oscilla¬ 
tions with large amplitudes. It can be seen from Fig. [3] that oscillations of the mode a(t)(f > 1 
with symmetry group Oh are not strictly sinusoidal. The frequency v of this mode depends 
on its amplitude A, and we depict the dependence v(A) in Fig. SI For small amplitudes the 
frequency changes slightly, while for large amplitudes the considerable decreasing of u(A) 
with increasing of the amplitude takes place. Therefore, a{t)(f ) 1 vibrations demonstrate soft 
type of nonlinearity. In the limit A —> 0, the frequency of the discussed mode tends approx¬ 


imately to 710 cm 1 .This value may be compared with the frequencies 718 cm 1 24] and 


723.5 cnr 1 25] for the mode zq[A 9 ], which were found with the aid of DM 0 I 3 code, and 
with the experimental value 772.27 cnr 1 . Possible reasons for discrepancy between above 
given values of z/i[A 9 ] will be analysed elsewhere. In this paper, we focus only on confir¬ 
mation of the general bush theory of nonlinear vibrations in physical systems with discrete 
symmetry using SFq molecule as a simple example. According to this theory in the case of 
large amplitudes, we should observe a very de fini te interaction between normal modes which 
were negligible for the small-amplitude vibrations. 

In Fig. 0 we see interaction between the vibrational modes b{t)<p 2 and a(t)(f ) 1 whose sym¬ 
metries are D^h and Oh , respectively. Indeed, at the initial instant t — 0 only mode b{t)(f > 2 
(root mode) was excited. However, the excitation from this mode is gradually transferred to 
the mode a(t)(f ) 1 which is the secondary mode [its symmetry is higher than that of the mode 
b(t)<p 2 \. The interaction between the above modes is very important since the amplitude of 
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the secondary mode reaches 30%, and even more, as compared with that of the root mode. 
Moreover, dynamics of the secondary mode is radically different from the simple sinusoidal 
oscillations, and this difference becomes more visible with increasing amplitude of the root 
mode. 

The main point is that excitation from the root mode is transferred to only one of all 
other modes, and exactly to the mode which must be excited in accordance with the bush 
theory (excitation can be transferred only to the modes with higher or equal symmetry). 

In our case, the secondary mode a(t)(f> 1 is associated with the two-dimensional irrep with 
conventional optics symbol E g . We have obtained the frequency v[D^ of this mode for 
small-amplitude vibrations with 9% error according to its experimental value, while the 
root-mode frequency u[Oh] was calculated with 8 % error. However, it is interesting to note 
that the ratio v[Oh] \ v[D±h\ found in our calculations is equal to 1 . 2 , and this value almost 
exactly coincides with that obtained from the experimental data. 

All above discussed manifestations of nonlinearity in dynamics of the two-dimensional 
bush B[D^h\ can be also found in time-evolution of the three-dimensional bush B[C\ V ] whose 
dynamics is shown in Fig. [HJ From this figure one can see that two secondary modes a{t)4> l 
and b(t)<fi 2 only these modes!) are involved into vibrational regime due to the interaction 
with the root mode c(f )</> 3 associated with the three-dimensional irrep iq ]. 

According to [25), the infrared absorption spectrum of SF 6 is composed of two bands 
whose maxima are located at v 3 = 938 cm -1 and zq = 567 cm -1 . Both of these normal 
modes are triple degenerate since they belong to the irrep F*. Our rough estimate of 
c(t)(f) 3 - mode frequency in the case of small amplitude vibrations gives u\C^ = 530 cm -1 
that allows us to identify this frequency with 1 / 4 . 

The Figs. [5] and 0 clearly show that the bush with dimension m > 1 represents a quasi- 
periodic dynamical object because it is a superposition of modes with different frequencies. 

In conclusion, let us note that only one mode, c(t)<p 3 , belonging to the three-dimensional 
irrep iq] contributes to the bush B[C\ V ] . What would happen if we initially excite another 
mode of the irrep iq) ? It is easy to understand that the result of this action will be emergence 
of another bush which turns out to be dynamically equivalent to the bush B[C^ V \. Such 
bushes we call “dynamical domains” of one and the same bush, in full accordance with the 
corresponding term of the theory of phase transitions in crystals. 

In the case of the irrep there exist three dynamical domains of the bush B[C^ V ] 
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which differ only in directions of the atomic vibrations. The axis of symmetry of the above 
examined bush B[C± V ] coincides with the coordinate axis Z , while the axes of symmetry of 
two other its domains coincide with X and Y coordinate axes. 


IV. CONCLUSION 


The main goal of this paper was to verify the validity of the group-theoretical results 
obtained by the theory of bushes of nonlinear normal modes in physical systems with discrete 
symmetry using as an example SF§ molecule. We have also investigated some properties of 
one-, two- and three-modes exact nonlinear vibrational regimes (bushes of NNMs) in this 
molecule with the aid of ab initial simulations based on the density function theory. For 
these calculations we have used ABINIT code 22], 


Some problems concerning the above bushes in SFq molecule turn out to be beyond the 
scope of the present paper. They will be studied in further works. Some of these problems 
are: 


1. Stability of the bushes of NNMs in SF 6 molecule which were obtained by the ABINIT 
code; 

2. The method for exciting these bushes by external fields; 

3. The possibility to describe large nonlinear vibrations in SF$ in the framework of 
mass-point models whose particles interact via pair phenomenological potentials; 

4. The possibility of applying the theory of bushes of nonlinear normal modes for studying 
some structural phase transitions whose nature relates to vibration of octahedron clusters 
in perovskite-like crystals. 
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FIGURES 



FIG. 1. Octahedral molecule. 
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b) 

FIG. 2. Displacement patterns for modes 4> l , <f> 2 and <f> 3 . 


22 




















FIG. 3. Vibrations of the bush B[Oh ] with different amplitudes. Note that the scale on the vertical 
axis is different for the cases a, b, c. In all figures, the atomic displacement are given in Bohrs, 
while time t is given in the atomic units ( a.u .). 
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FIG. 4. Dependence of the frequency v of the breathing mode on its amplitude A. 
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FIG. 5. Dynamics of the bush B[D^ with different amplitudes of the root mode b(t). 
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max|a(t)| 



FIG. 6. Dependence of max|ai(t)| on max 1 0 , 2 (t)|- 



FIG. 7. Approximate dynamics of the bush B[D^\. 
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FIG. 8. Vibrations of the bush B[C± V \ with different amplitudes of the root mode c(t). 
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